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It is widely believed that chiral symmetry is spontaneously broken at zero temperature in the 
strong coupling limit of staggered fermions, for any number of colors and flavors. Using Monte Carlo 
simulations, we show that this conventional wisdom, based on a mean-held analysis, is wrong. For 
sufficiently many fundamental flavors, chiral symmetry is restored via a bulk, first-order transition. 
This chirally symmetric phase appears to be analytically connected with the expected conformal 
window of many-flavor continuum QCD. We perform simulations in the chirally symmetric phase 
at zero quark mass for various system sizes L, and measure the torelon mass, the Dirac spectrum 
and the hadron spectrum. All masses go to zero with 1/L. L is hence the only infrared length 
scale. Thus, the strong-coupling chirally restored phase appears as a convenient laboratory to study 
IR-conformality. Finally, we present a conjecture for the phase diagram of lattice QCD as a function 
of the bare coupling and the number of quark flavors. 

PACS numbers: 11.30.Rd, 12.38.Gc, 12.60.-i 



The possibility that the Higgs boson could be a composite bound-state in a high-energy Technicolor theory [TJ 
has generated considerable interest, especially in the lattice community. In particular, the requirement that the 
Technicolor theory be "walking" 2J, in order to accommodate stringent bounds on flavor-changing neutral currents, 
has been the driving motivation behind several large-scale computer simulation efforts to determine the possible 
combinations of gauge groups and fermion contents leading to a conformal window. 

To determine via lattice Monte Carlo simulations whether a given theory is inside the conformal window is partic- 
ularly challenging, because it involves a triple difficulty: in order to identify (or not) an infrared fixed point (IRFP) 
which is the signature of a theory inside the conformal window, one must probe the extreme infrared properties of the 
theory, while at the same time taking the continuum limit of the lattice discretization, and controlling the limit when 
the quarks become massless. This compounded difficulty may explain why, in spite of considerable efforts, there is no 
consensus yet on the minimum number Nf* of quark flavors needed for QCD to be inside the conformal window [3]. 
A numerical demonstration of walking has been provided only recently, in a toy model, the 2-d 0(3) model at vacuum 
angle 9 sa tt [3] . 

Here, we relax the demand that results should be obtained in the continuum limit. On a coarse lattice, long-distance 
properties can be studied more economically. While such properties may differ from those of the corresponding 
continuum theory, it may still be instructive to consider the possible existence of an IRFP for a discretized lattice 
theory. The phase diagram of SU(A r c ) gauge theory with Nf fundamental fermions, as a function of Nf and the 
bare gauge coupling, has been predicted in the celebrated Ref. [5] , which serves as a guide to understand the results 
of Monte Carlo simulations performed at finite bare coupling. It is important to confront these predictions with 
uncontroversial numerical evidence. Therefore, we start our investigation by considering the strong coupling limit, 
where the lattice is maximally coarse. 

Note that we consider standard staggered fermions, and (away from the strong-coupling limit) the standard plaquette 
action. Other discretizations could lead to different results, since only the continuum limit is universal. 

The conventional wisdom for strong coupling QCD with staggered fermions is that chiral symmetry remains always 
broken at zero temperature, regardless of the number of colors and flavors. This belief is based on mean-field analyses 
performed in some of the earliest papers on lattice QCD. In particular, it was shown in [6] that at leading order in 
a l/d expansion, the chiral condensate has a value independent of Nf and N c , but depends only on the number d of 
spatial dimensions: 



Chiral symmetry may be restored by increasing the temperature T. Following the approach of [7], where explicit 
results are provided for a few small values of N c and Nf, we calculated the chiral restoration temperature aT c and 
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found that it is indeed non-zero for all Nf , and independent of N c to leading order in 1/ Nf : 
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Hence chiral symmetry will never be restored at zero temperature, according to the mean-field analysis. Since mean- 
field theory is expected to work well when the number of d.o.f. per site is large (e.g. providing exact results in the 
Gross-Neveu model for Nf —> oo), there was no reason to doubt the validity of this finding. Besides, it is in accord 
with the intuition that the gauge field is maximally disordered in the strong-coupling limit, and that this disorder will 
drive chiral symmetry breaking. 

On the other hand, one naively should expect the above disorder to be modified by dynamical fcrmions, which have 
an ordering effect. Indeed, the loop expansion of the determinant shows that the fermionic effective action induced by 
dynamical fermions, S'cff — — log det(0 + m q ), starts with a positive plaquette coupling A/3, proportional to 1/m^ for 
heavy quarks, which has been studied numerically in [8. . Clearly, for Nf flavors the effective action is proportional 
to Nf, and A/3 grows proportionally. This plaquette term suppresses fluctuations in the gauge field, which suggests 
that chiral symmetry restoration might take place for sufficiently large Nf. 



II. MONTE CARLO RESULTS 
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FIG. 1: The chiral condensate at strong coupling, j3 = 0, in the (Nf,am q ) plane, for 4 4 (left) and 6 4 (right) lattices. 
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FIG. 2: The chiral condensate at weaker coupling, /3 = 5, in the (Nf,am q ) plane, for 4 4 (left) and 6 4 (right) lattices. 



The only way to resolve this puzzle is to perform Monte Carlo simulations in the strong coupling limit of staggered 
fermions, to detect a possible chiral symmetry restoration for sufficiently large Nf. These simulations are straight- 
forward, using the standard Hybrid Monte Carlo algorithm. As expected, the effect of increasing Nf on the chiral 
condensate is to reduce its magnitude. But it came as a surprise to find that the chiral condensate vanishes via a strong 
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FIG. 3: The chiral condensate at am q = 0.025 in the (Nf,f3) plane, for 4 4 (left) and 6 4 (right) lattices. The phase transition 
remains strongly first order at weaker coupling. The contour plot indicates the qualitative behaviour of the phase boundary, 
extending to weaker coupling. 

first-order transition at Nf C ~ 52 continuum flavors in the chiral limit (i.e. Nf C /A ~ 13 staggered fermion fields). 
In the broken phase, the chiral condensate remains almost constant. It vanishes in the chiral limit due to finite-size 
effects only. In contrast, in the chirally restored phase the condensate is caused by explicit symmetry breaking and is 
proportional to the quark mass. This is illustrated Fig. [TJ where the condensate is shown as a function of Nf and bare 
quark mass {am q ). Moreover, this Nf -driven transition turns out to be a bulk, zero-temperature transition, which 
can be seen by the fact that finite-size effects on the phase boundary are small when comparing two different system 
sizes 4 4 and 6 4 , as shown Fig. [l] (left and right). 

One also observes in these figures that the critical number of flavors increases with the quark mass. This is easy 
to understand: heavier quarks have a weaker ordering effect, so that the induced plaquette coupling A/3 decreases if 
one keeps Nf fixed. It takes more flavors to keep the system chirally symmetric. Hence, Nf C increases, and for heavy 
quarks should obey Nf C cx (am q ) 4 . 

One can now go back to the mean-field treatment and trace the origin of its failure. Two kinds of terms at least are 
neglected: (i) multiple meson hopping along a given link, and (ii) baryon loops. These terms amount to corrections 
0(Nf/N c ) and O(NffcP), respectively, where Nf = Nf /A is the number of staggered fields, and is normally set to I 
in the mean-field treatment. Here, we consider Nf > 13, and the previously neglected corrections become dominant. 
The conventional wisdom that chiral symmetry is always broken at T = in the strong-coupling limit comes from 
mistakenly applying the lowest-order mean-field approximation in a regime where it is invalid. 

Having established an iV/-driven phase transition in the strong-coupling limit, we may consider its impact on the 
lattice theory at non-zero lattice gauge coupling /3 as well. Since the transition is strongly first order, it has to persist 
for some range in /3 at least. Hence we have compared the strong coupling phase diagram with the phase diagram at 
weaker coupling ft — 5, illustrated Fig. [5] We find a similar qualitative behavior, but with Nf C drastically reduced to 
O(10). Finite-size effects are more pronounced, but the transition still seems to be a first-order bulk transition. 

In fact, we find a smooth variation of the Nf -driven transition with /3 at a given small quark mass am q = 0.025, 
as shown Fig. [3] The transition extends to weak coupling, at least to (3 = 5, and remains strongly first-order. Thus, 
it is plausible that this transition, which separates a chirally broken (small Nf) and a chirally symmetric (large Nf) 
phase, persists all the way to the j3 — ► oo continuum limit, where it is to be identified with the transition at Nf — Nf* 
between the chirally broken and the IR-conformal, chirally symmetric phase. In other words, our chirally restored 
phase may be analytically connected to the conformal window in the continuum limit, because we do not observe any 
additional non-analyticity as /3 is increased. 

This possibility motivates our study of the properties of the strong-coupling chirally symmetric phase, looking for 
tests of IR-conformality. 

III. LOOKING FOR CONFORMALITY IN THE CHIRALLY SYMMETRIC PHASE 

It is natural to ask whether the chirally restored phase is connected to the conformal window, i.e. whether the 
chirally restored phase at strong coupling is also IR-conformal. And if this is indeed the case, the next obvious 
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question is whether this IR-conformal phase is trivial, ie. whether the IR fixed point coupling is zero or not. In 
this section we present measurements of gluonic and fermionic observables chosen to address these questions: the 
torelon mass from which we define a running coupling, the Dirac eigenvalue spectrum, and the hadron spectrum. Our 
results support the following conclusion: the strong-coupling chirally symmetric phase is indeed IR-conformal, and it 
is non-trivial. 

The simulations performed here are all in the chirally symmetric phase at zero plaquette coupling, with Nf = 56 
and 96 continuum flavors, and with lattices of size 4 3 x 16, 6 3 x 16, 8 3 x 16, 10 3 x 20 and 12 3 x 24. The quark mass 
is set exactly to zero unless specified otherwise. We will see below that the Dirac operator has a spectral gap in the 
symmetric phase, which makes a study of the massless theory within reach of modest computer resources. Moreover, 
having one infrared scale, the system size L, rather than two scales (L and l/m q ) is of great advantage when analyzing 
the results. 

Let us mention the average plaquette values which we measure: ~ 0.35 and ss 0.52 for Nf — 56 and 96, respectively 
(normalized to 1 for the free field). So we are very far from a plaquette value of 0, corresponding to maximally 
disordered gauge fields and achieved for Nf — 0: the ordering effect of the dynamical fermions plays a dominant role 
in our case, and the vanishing of the plaquette coupling is not associated with special properties. 



A. Characterizing the chirally restored phase: (I) The Torelon Mass 
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FIG. 4: Torelon correlator 

The "torelon" is a gluonic excitation which is topologically non-trivial: it is excited by any Wilson loop which wraps 
around the spatial boundary in one direction, for instance, as illustrated Fig. [4j 

L-l 

T i (t)= r Trl[U i (x + ki,t), (3) 

where i = 1,2,3 is one of the spatial directions and i the unit vector in that direction. 1 We extract the mass of this 
excitation from the exponential decay of the correlator (Xi(O)Tj(t)). To suppress excited states, we smear the links 
within each time-slice before constructing Ti. This observable has been used for a long time to extract the string 
tension a in Yang-Mills theories [9j: it can be viewed as a loop of gluonic string, whose energy rnriL) grows with 
its length as ah. So our initial, naive expectation was to measure a mass m<r(L) growing with L, until perhaps the 
string would break due to fermion-pair creation. 

This is not at all what we observed. The dimensionless quantity which we measure, amxiL), decreases on larger 
lattices corresponding to a larger ratio L/a. Clearly, our theory is not confining. Moreover, as shown Fig. [5] the 
combination LrriT(L) is approximately constant as L is increased. So the torelon mass varies as 1/L (actually, for 
small L it initially decreases even faster with L as seen in the figure). Thus, there is no intrinsic mass scale which 
appears in this channel: the torelon mass is set by the system size L. This remarkable result is our first evidence that 
our theory is IR-conformal. 



1 We initially adopted periodic boundary conditions in space and anti-periodic in time for the fermion fields. However, the dynamical 
fermions drive T; to negative values in this case, with a Z2 degeneracy between the two complex Zz sectors. We occasionally observed 
tunneling of Ti between these two sectors, and long metastabilities. For this reason, we changed the spatial boundary conditions to 
anti-periodic, which makes (Ti) real positive. 



5 



(Torelon mass x L) vs L, N f =56, m q =0 
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FIG. 5: The torelon mass rriT(L) multiplied by L versus L, left: Nf = 56, right: Nf = 96 (for anti-periodic boundary 
conditions). 
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FIG. 6: The running coupling g(L) defined via the temperature-dependence of the Debye mass, which is identified with the 
torelon mass, for Nf = 56 and 96. For each Nf, the curve is a fit to a constant plus (a/L) 2 corrections of the 4 largest volumes. 
The larger Nf has a smaller coupling. 



We actually combined the Ti(t),i £ {1,2,3} into two representations of the cubic discrete rotation group: the 
(corresponding to a + representation of 0(3)) and the E + (corresponding to a 2 + representation of 0(3)). The mass 
of the E + seems to be slightly smaller, as observed in small- volume analytic Yang-Mills calculations 10J. 

Now, by relabelling the spatial direction i as the imaginary time direction, one realizes that we are measuring the 
correlation of two time-like Polyakov loops, whose decay rate is governed by the Debye mass, given perturbatively at 
lowest order by 



»'vlT) = -2 !l T\l^ + %- (4) 



This expression allows us to define a running coupling g(L) via 



9(L) = (5) 

3 ' 6 

and we see that, in this scheme, our running coupling seems to go to a non-zero constant as L increases (although 
one cannot exclude, of course, that it slowly goes to zero). Therefore, we have numerical evidence supporting the 
view that our strong-coupling, chirally symmetric theory is IR-conformal and non-trivial. 
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Interestingly, the extracted value of g(L) approaches ~ 0.95 and ~ 0.80 for Nf — 56 and 96 respectively. So the 
IR fixed-point coupling value decreases as Nf increases. This is what one would expect: as Nf keeps increasing, the 
ordering effect of the fermions increases, and all Wilson loops are driven towards I, their free field value. At the same 
time, any definition of a running coupling will approach zero. The theory becomes trivial for Nf — > oo, even in the 
strong-coupling limit. We will come back to this point in Sec. IV. 



Integrated eigenvalue density, 10 , m =0 



Integrated eigenvalue density, rescaling the infrared, m q =0, N f =56 
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FIG. 7: The integrated eigenvalue density. Left: comparison of Nf = (quenched) with Nf = 56 in the chirally restored phase, 
where only the latter shows a spectral gap. Right: The rescaled spectral gap, indicating 1/L scaling. 




FIG. 8: The anomalous dimension from fits to the spectral gap for Nf = 56, 96 with periodic and anti-periodic boundary 
conditions. Fitted exponents to the three largest volumes are 7* ~ 0.26 for Nf = 56 and 7* ~ 0.38 for Nf = 96. 



B. Characterizing the chirally restored phase: (II) Dirac Spectrum 



We now turn to fermionic properties, starting with the spectrum of the Dirac operator. We have analyzed the 
Dirac eigenvalue spectrum of the configurations in our Monte Carlo ensembles, using a Lanczos algorithm to obtain 
an approximation of the whole spectral density and an Arnoldi method to extract the smallest eigenvalues to high 
accuracy. The observable shown in Fig. [7] is the integrated eigenvalue density, defined as: 

This function of A counts the fraction of eigenvalues smaller than A. Its derivative is simply the eigenvalue density 
p(X). We first compare this observable on quenched configurations (Nf — 0) and in the chirally symmetric phase 
(N f = 56). In Fig. W\ (left), 10 curves for 10 configurations are superimposed: variations in the spectrum are very small. 
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FIG. 9: Left: meson masses and m p for various system sizes L, showing a fast decrease of the masses as L increases, while 
the mass ratio remains more or less constant. Right: pion mass as a function of with the three largest volumes fitted by 
the ansatz (1/L) ^ 1+7 \ yielding 7* ~ 1.0 and 0.4 for N; = 56 and 96, respectively. 

We observe that Nf = and Nf — 56 spectra are similar in the ultraviolet, but differ in the infrared, as illustrated in 
the inset. The Nf — curve starts linearly from the origin, reflecting an eigenvalue density approximately constant 
near A = 0. On the contrary, the integrated eigenvalue density for Nf — 56 shows a spectral gap for small eigenvalues, 
which is of course consistent with chiral symmetry restoration according to the Banks-Casher relation, since p(Q) = 0. 

The crucial question is on which scale does this spectral gap depend. To answer this question, in Fig. [7] (right) we 
compare the integrated eigenvalue density at Nt = 56 for various lattice volumes, L = 4,6, 8, 10 and 12. As evidenced 
in the inset, we find that the spectral gap scales oc 1/L to a good approximation, which is a strong indication that 
our theory is IR-conformal 2 : There does not seem to be any length scale in the chirally restored phase other than the 
box size L. 

Actually, small deviations from 1/L scaling allow us, in principle, to extract the anomalous mass dimension 7*. 
We make such an attempt in Fig. [8j where the gap has been multiplied by L already: deviations from a constant are 
indicative of anomalous dimension, provided other corrections 0((a/ L) 2 ) are negligible. The effect of a finite system 
size L on the Dirac spectrum has not been analyzed yet. We have simply considered that the infrared conformal 
symmetry is explicitly broken by the infrared scale 1/L, which is the analogue of an explicit breaking by a quark mass 
m q . Consequently, we expect the mass gap to behave as {1/ L) x ^ 1+1 '\ A crude fit based on our 3 largest volumes 
gives 7* ~ 0.26 and 0.38 for Nf = 56 and 96, respectively. 

C. Characterizing the chirally restored phase: (III) Hadron Masses 

Finally, we turn to hadron masses measured on our Nf = 56 and Nf = 96 Monte Carlo ensembles. Even though 
the quark mass is zero in these ensembles, we observe non-zero hadron masses. As expected, parity partners are 
degenerate since chiral symmetry is restored. Now, if our theory is IR-conformal, the masses which we measure are 
exclusively due to finite-size effects: all masses should go to zero as the lattice size L is increased. This is what we 
observe, as shown in Fig. [9] left: the "pion" and "rho" masses both decrease by a factor ~ 2 as the lattice size is 
increased from L — 4 to 12. Notice however that the smallest mass is still > 0.6, which is not very light. 

Furthermore, one generally expects that the approach to zero should be the same for all hadrons, so that mass ratios 
should remain constant as L —> 00. Note that there may be exceptions to this "rule": in the 2d 0(3) sigma model 
near 9 = ir, the mass of the 0(3) singlet state approaches zero faster than that of the 0(3) triplet as 9 approaches 
7T [11] . Here, our Fig. [9] left would show all data points aligned on "rays" going through the origin if mass ratios were 
constant. One can see deviations from this behavior, which perhaps are caused by the not-too-light masses which 
we measure. Another possible cause is technical: as L is increased, the groundstate masses in each channel decrease, 



2 If one would take the limits L — > 00 first, then m q — > 0, the expected spectral density for a conformal theory would be p(A) ~ 
\(3—~t )/(i+7 ) Here, we take the opposite order of limits. 



but so do also the mass differences between groundstate and excited states. It becomes more difficult to extract the 
groundstate mass, and our lattices are likely too short to properly control this source of systematic error. 

Nevertheless, we show in Fig. [9] right the mass of the "pion", in which we have the most confidence, as a function 
of 1/L. Since l/L breaks the conformal symmetry much like a quark mass m q would, we expect that the pion mass 
should scale the same way, namely as (1/L) 1 /( 1+7 \ if L is large enough. A fit to our three largest system sizes gives 
7* ~ 1.0 and 0.4 for Nf = 56 and 96, respectively. 



IV. THE CONJECTURED PHASE DIAGRAM 



Using both gluonic and fermionic observables, we have presented evidence that the chirally restored phase at 
strong coupling is IR-conformal and non-trivial, and speculated about a connection to the conformal window in the 
continuum. We now want to propose a phase diagram sketched in Fig. |10[ as a function of the plaquette coupling 
/3 = 6/(?o an d of the number of fundamental flavors Nf. This conjectured phase diagram can be compared with, 
e.g., those of Ref. 0[T3]: one can see substantial differences. The essential feature of our phase diagram is that the 
P = IR-conformal phase is analytically connected with the weak-coupling, continuum IR-conformal phase. This is 
the simplest scenario, supported by our exploratory scan in j3 shown Fig. A single phase transition line Nf c (j3) 
separates the region of broken chiral symmetry at small Nf from the chirally symmetric region at large Nf. The 
transition is first-order, at least for some range of /3 starting from zero. 

Furthermore, based on our results for the running coupling Fig. [6] we conjecture that for any finite j3 and Nf > 
Nf c (/3), large- Nf lattice QCD is IR-conformal, with a non-trivial fixed-point coupling g* > 0. This value changes 
continuously with (f3,Nf), reaching the value zero for j3 = oo, Nt > 33/2 and for Nf — oo V/3, as indicated Fig. 10 
left by a thick dotted line, g* grows as one moves away from this dotted boundary towards the phase transition line. 

An alternative scenario would be that the running coupling in Fig. [6] slowly approaches zero, and g* = for Nf > 
33/2 for any j3 in the chirally symmetric phase. If the basin of attraction of the weak coupling trivial fixed point would 
extend all the way to the strong coupling limit, one should observe for the running coupling g 2 (L) ~ l/log(L/Lo), 
with L ~ O(a). This possibility is indeed consistent with our numerical results Fig. [6] for the running coupling, but 
should give an anomalous mass dimension 7* = 0, which is not favored by our other measurements. More careful, 
large-scale simulations are necessary to settle this issue. 

Finally, one may consider the line g* = go, with go = (2iV c //3) 1/ ' 2 , where the IR fixed-point coupling has the same 
value as the bare coupling, so that the coupling does not run as a function of the renormalization scale. This line 
starts at the point (/3 — 00, Nf = 33/2) where g* = go = 0. Its precise location depends on the chosen renormalization 
scheme. It is not associated with any kind of singularity of the free energy. There is no phase transition along this 
line. 

We have determined the phase diagram in the strong-coupling region only. Studying the continuum limit is of 
course much more difficult, due to the large lattices that have to be used in order to control the finite size effects 
and the difficult control of lattice artifacts. We would like to suggest that the strong-coupling limit may represent an 
advantageous "poor man's laboratory" for the study of 4d IR-conformal gauge theories. In particular, as illustrated 
Fig.[l0]right, the range of scales over which conformal invariance applies for a given computing effort is greatly reduced 
at weak coupling: there, for a given lattice size N 4 , the scales are ordered as follows: 

a<Cl/A«£ = Aa. (7) 

where A is the scale generated by the asymptotically free gauge dynamics. In contrast, at strong coupling, where the 
lattice becomes maximally coarse, the hierarchy is: 

a- 1/A< L = Na. (8) 

Hence the dynamical range of conformal invariance, characterized by the product LA, is maximized at (3 = for a 
given lattice size N = L/a. 



V. CONCLUSION 



We have shown that for j3 = 0, contrary to common wisdom, there exists a strongly first-order, iV/-driven bulk 
transition to a chirally symmetric phase. In the chiral limit, the transition occurs for Nf c = 52(4) continuum flavors. 
This finding is in contrast to the mean-field prediction, whose failure can be traced back to approximations relying on 
Nf being small. Clearly, the conventional, automatic association of the strong-coupling limit with confinement and 
chiral symmetry breaking is too naive. Furthermore, the chirally restored phase extends to weak coupling. 
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FIG. 10: Left: conjectured phase diagram in the (/3,Nf) plane. A single phase transition line separates the chirally broken 
phase from the chirally symmetric, IR-conformal phase. The thick dotted line indicates trivial theories. Right: the ordering of 
scales at weak coupling (top) and strong coupling (bottom), indicating that the range of conformal invariance is larger in the 
latter case. 



We have also shown numerical evidence that the (3 — chirally restored phase of "large- Nf QCD" is IR-conformal, 
with a non-trivial, iVf -dependent value of the IR fixed-point coupling. We conclude that the strong-coupling limit is 
the laboratory of choice to study a 4d IR-conformal gauge theory. Simulations at large Nf and zero quark mass can 
be performed without too much computational effort since a gap appears in the Dirac spectrum. As Nf increases, 
the spectral gap increases, the average plaquette approaches 1, and the fixed-point coupling approaches 0. Setting 
the quark mass to zero eliminates one IR scale, leaving the system size L as the only remaining one. This greatly 
simplifies the analysis of simulation results. 

Since we have not observed any evidence for an additional T = phase transition as /3 is increased, we speculate 
that the strong coupling chirally symmetric, IR-conformal phase is analytically connected with the continuum 
IR-conformal phase. 



One may ask how robust these statements are with respect to the particular discretization of the Dirac operator and 
the gauge action. While the quantitative details of the phase transition Nf C (/3) will surely change, we think that the 
qualitative features will remain. Chiral symmetry breaking at strong coupling, for small Nf, is a general consequence 
of the disorder in the gauge field. The ordering effect of fermions also is generic. So we do expect a bulk transition, 
gcncrically of first-order, as a function of Nf in the strong-coupling limit. Actually, such a transition was observed for 
Wilson fermions in Ref. [T3]- At intermediate coupling, additional transitions may appear depending on the lattice 
action. Interestingly, a first-order transition to a chirally broken phase as /3 is reduced has been observed many times, 
for various lattice actions [141 - 1X8] . These transitions were observed for some fixed value of Nf. Here, we simply put 
all these earlier observations together. It is interesting that this phase transition is consistently of first-order. If the 
first-order nature persists all the way to the continuum limit, then walking dynamics will not be observed, and the 
transition to the conformal window will be characterized by "jumping dynamics", as proposed by Sannino |19j . 

One may wonder what happens to the (/3, Nf) phase diagram as the gauge group or the fermion content is changed. 
For 4d compact U(l), the change would not be dramatic, because the strong-coupling behavior is much the same as 
for SU(3); our first-order transition line would end at (j3 ss 1.01, = 0) on the horizontal axis rather than on the 
vertical axis, and the region of triviality would cover the whole chirally symmetric phase, except for Nf — 0. For 
SU(2) or SU(3) with adjoint fermions, the change would be more significant: in the strong-coupling limit, increasing 
Nf would order the plaquette in the adjoint representation, not in the fundamental. Center monopoles would likely 
condense |20j . and might delay or prevent the restoration of chiral symmetry. 

Finally, there are many directions in which to extend this exploratory study. More observables, like the static 
potential and the Fredenhagen-Marcu order parameter, should be studied. Also, and to make contact with other 
numerical studies, a mass deformation could be introduced. As a first step in this direction, we show in Fig. [TT] 
the quark mass dependence of the chiral condensate. This figure shows all the technical difficulties associated with 
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FIG. 11: Left: The chiral condensate as a function of the quark mass, showing clearly a linear relation as appropriate for the 
chirally symmetric phase. Deviations from linearity can be used to determine the mass anomalous dimension according to 
Eq.[9] Right: the same data, now divided by (am q ) to emphasize deviations from linear behaviour. The fitted values of 7* 
depend on the fitting range (delimited by the vertical line), but tend to be negative (e.g. 7* ~ —0.5 for L — 10). 



extracting the anomalous dimension 7*: Heavier fermions have less ordering effect, which triggers a phase transition 
back into the chirally broken phase for some critical fermion mass. Finite-size effects associated with that transition 
should not be included in the determination of 7*. Moreover, the non-anomalous contributions (ipip) — c±m q + Civn? q 

easily overwhelm the anomalous 7 • ) ^ 1+7 Actually, fits based on the ansatz 

3-7* 

= c\m q + c 2 m q + ~ l " (9) 

favor negative values for 7*, which crucially depend on the fitting range and the included analytic contributions. It 
seems that analyzing the Dirac spectrum provides better control of the systematics, as emphasized in |21j . 
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